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Abstract 

Recently, Marzlin and Sanders (2004) demonstrated an inconsis- 
tency when the adiabatic approximation was applied to specific, "in- 
verse" time-evolving systems. Following that, Tong et al. (2005) 
showed that the widely used traditional adiabatic conditions are in- 
sufficient to guarantee the validity of the adiabatic approximation for 
this class of systems. In this article we explore the origin of these obser- 
vations by a perturbative approach and find that in first order approxi- 
mation certain nonzero terms appear in the solution which gives rise to 
the breakdown of the adiabatic approximation (despite the fact that 
the traditional adiabatic conditions are satisfied). We argue that in 
this case the Hamiltonian of Marzlin and Sanders' inverse time evolv- 
ing system cannot be written in terms of t/T, where T denotes the 
total evolution time. It is further demonstrated that the new qualita- 
tive adiabatic condition of Ye et al. (2005) performs well in some cases 
when the traditional conditions fail to describe properly non-adiabatic 
evolution. 



1 Introduction 

The quantum adiabatic theorem is one of the oldest fundamental results in 
quantum physics [1] . It concerns development of systems where the nonde- 
generate Hamiltonian evolves slowly in time. In the limit when the change 
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of Hamiltonian H(t) is made infinitely slow the system, which started from 
one of the eigenstates of H(0), passes through the corresponding instanta- 
neous eigenstate of H(t) [2, 3, 4, 5]. The adiabatic theorem underlies the 
adiabatic approximation scheme, which states that if the Hamiltonian H(t) 
evolves slowly enough by satisfying the adiabatic condition in time interval 
t e [0, T], then the evolving state of the system will remain close to its 
instantaneous eigenstate up to a multiplicative phase factor in the interval 
[0, T\. The adiabatic approximation has potential applications in several 
areas of physics such as the Landau-Zener transition in molecular physics 
[6], Gell-Mann-Low theorem in quantum field theory [7] or in the lore of 
Berry phase [8]. Recently, with the emergence of the new field of quantum 
information theory , new interest has arisen in the application of the quan- 
tum adiabatic theorem [9] . An alternative scheme appeared beside the usual 
quantum algorithms, based on an adiabatic time evolution of the state of the 
Hamiltonian where the final Hamiltonian encodes the solution for a given 
problem [10]. In another application, the Berry phase has been proposed 
to perform quantum information processing tasks [11]. Both of these con- 
cepts exploit the circumstance, that the adiabatically evolving ground state 
is very robust against decoherence and small perturbations [12]. It is thus 
important to explore the limits of the adiabatic conditions in order to un- 
derstand better when the evolution of a quantum system can be considered 
adiabatic. 

In a recent Letter Marzlin and Sanders [13] demonstrated on a dual pair 
of systems (called a-system and b-system), that if the evolution operator 
of the b-system is the Hermitian conjugate of an adiabatically evolving a- 
system, then the application of the adiabatic approximation in the b-system 
can lead to contradiction. Tong et al. [14] explained this inconsistency by 
pointing out that the widely used traditional adiabatic conditions are not 
sufficient to guarantee adiabaticity. 

In the present paper our aim is to find the common root of the problems 
by analyzing the time evolution of the dual systems in a perturbative man- 
ner. In a first order approximation we obtain certain non-vanishing terms 
in the regime where the traditional adiabatic conditions hold. We point out 
that these nonzero terms are responsible for the violation of the adiabatic 
approximation. We further argue that in this case the Hamiltonian of the 
b-system cannot be written in terms of t/T, where T denotes the total evo- 
lution time. In addition, an explicit example is provided which on one hand 
illustrates our arguments, and on the other hand allows a new adiabatic 
condition of Ye et al. [16] to be tested on. 
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2 Analysis of dual pair of systems 

2.1 Expansion of state 

Consider a closed iV-dimensional quantum system in a state \4>{t)), which 
evolves through the time-dependent Schrodinger equation 

where H (t) denotes the time-dependent non-degenerate Hamiltonian of the 
system and we set h = 1. Let us introduce the normalized time s by the 
variable-transformation t = sT, < s < 1, and rewrite equation (1) as 

i^(s,T)) =TH(s,T)^(s,T)) , (2) 

where T denotes the total evolution time. The instantaneous eigenstates 
\E n (s,T)} of the Hamiltonian H(s,T) satisfy 

H(s,T)\E n (s,T)) = E n (s,T)\E n (s,T)) , n = l,...,N, (3) 

and the elements of the non-adiabatic coupling matrix are defined by 

r nk (s,T) = (E k (s,T)\-^-E n (s,T)) . (4) 

We now expand the state \ijj(s,T)) in the basis of the instantaneous eigen- 
states of H(s, T), 

N 

4>(s,T) = ]T Ms,T)e- tT Jo E "W ds '\E n (s,T)) . (5) 

n=l 

(We shall frequently neglect in the text, but not in the formulae, the variables 
(s,T) from the expressions following from equations (2,3,4), and the possible 
(s,T)-dependence is understood without denoting it.) 

2.2 Parallel transport 

Let us now introduce the following local phase change for the n'th instan- 
taneous eigenstate 

\E n (s,T)) = e i@ ^\E n (s,T)) , n = l,...,N, (6) 
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where @ n are real, (s,T)-dependent parameters. Plugging (6) into the defi- 
nition (4) the transformation formula for the non-adiabatic coupling terms 
[17] reads 

f nk (s, T) = e^^- Q ^r nk (s, T) + i^^-S nk . (7) 

This simple relation says that the diagonal element f nn is boosted with 
respect to T nn by idQ n /ds, while the non-diagonal elements of f take up 
a phase with respect to those of r. If we choose the phase Q n (s,T) = 
i Jq S T nn (s', T)ds' , then under relation (7) f nn becomes zero, and by definition 
(4) the nth eigenstate satisfies the parallel transport law (E n (s, T)\d/ds\E n (s, T)) 
0. Let us denote in this gauge the matrix elements of r by rj fe , then we get 

T ) = ei ° Tkk{s '' T) - Tnn {s '' T)ds 'r nk (s,T) , (8) 

for n 7^ k. Substituting \ifj(s,T)) expressed in the rotating frame (5) into 
(2), and performing some algebra we obtain for the complex amplitudes (f) n 
the following differential equation: 

^ = - E Ms,T)rl{s,T)e- iT r ^,T W (g) 

where g n k(s, T) = E n (s, T) — E k (s, T). Now let us define the matrix elements 
A nk (s,T) by 

A nk (s,T) = Tnk ^l\ . (10) 
9nk{s,T) 

The widely used, traditional condition for the adiabatic approximation, if 
the system starts its evolution in the n'th instantaneous eigenstate of the 
Hamiltonian H(t), is then encoded in the statement [3, 13, 14, 18, 19] 



{E k (t)\E n (t)) 



E k {t)-E n (t) 



<1, k^n, t€[0,T\, (11) 



which under the variable transformation t = sT is equivalent to the condition 
[15] 

\A nk (s,T)\ «T, k^n, s €[0,1]. (12) 
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2.3 a and b systems 



Let us now turn our attention to a pair of A-dimensional quantum systems 
[13, 14, 20], where the a-system defined by Hamiltonian H a (t) and its dual 
b-system with Hamiltonian H b (t) are related through the following formula: 

H\t) = -U a \t)H a {t)U a {t) , (13) 

where we assumed that the spectrum of H a (t) is entirely discrete and non- 
degenerate. Formula (13) implies U b (t) = U a ^(t) between the evolution 
operators of the a-system and the b-system, and links the eigenvalues and 
the eigenvectors of the dual systems through the equations E^(s,T) = 
-E%(s,T) (and equivalent^ g b nk (s,T) = -< fc (s,T)) and \E b n (s,T)) = (s,T)\E a n {s,T)) . 
According to Duki et al. [20] after some manipulation of these formulae one 
may arrive at the following correspondence between the matrix elements of 
r" for the dual pair of systems: 

r£( S ,T) = rS( S ,T)e-^/o S ^(^'. (14) 

It is noted that since according to the definition (13) the a-system and the 
b-system are interchangeable, the equations numbered (13 — 20)have dual 
pairs by exchanging the labels (a) and (b). Now let us write out equation 
(9) specifically for the b-system, 

MM = - VJ ^( 5 ,T)r£( 5 ,T)e--«*t.',^ , (15) 

nj^k 

By use of (14) this results in a simpler expression for the evolution equation 
of the b-system [20], 

^^ = -E^r)r; (s ,T). (16) 

nj^k 



2.4 Perturbational solution 

In the following we consider the above pair of equations (15,16) and solve 
each of them perturbatively. Let us consider that the initial state of the 
b-system is |£*(0,T)) with <$.(0, T) = 5 n k- As zeroth order solution, the 
amplitudes <p b k are assumed not to be evolving, i.e. (j) b k {s,T) <p b k (0,T), 
which constitutes the adiabatic approximation, and in turn by insertion 
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into the rhs of (15) and (16) yields for 4> k the parallel set of first order 
solutions shown below. Performing an integration by parts as well, the first 
order solution for (15) is given by 

4(s,T) = 4(0,T) + ^J2( P nk(s,T) + Q b nk (s,T)) , (17) 

n^k 

where the terms P b k and Q b nk are specified by 

P b k (s,T) = A b nk (0,T) - A b nk ( S ,T)e- lT Io^ s '^ ds ' , (18) 

and 

QUs.T) = /;e--/:'i<-".^-"^k^d,' , (19) 
for n / k. On the other hand, the first order solution of (16) is 

4(s,T) = #(0,T) + J2 f^ll (s',T)ds' . (20) 

(The first term on the rhs of (17) and of (20) corresponds to the zeroth 
order solution.) We note that equation (17) would be formally identical to 
Eq. (21) in ref. [15] if Sarandy et al. fixed the gauge to obey the parallel- 
transport law r nn = 0. However, the main difference lies in the fact that 
Sarandy et al. do not allow H b (s,T) to depend on the total evolution time 
T. By contrast, in our treatment the Hamiltonian of the b-system and all 
the other quantities derived from it accommodate a possible T-dependence. 
(It is noted that from the notion of T-dependence we exclude the trivial 
situation when T appears as T/T', where T' is different from but the same 
order as T, since here the presence of T is immaterial.) 

Now let us consider the case [15, 19] when the Hamiltonian of the a- 
system H a is a (continuous) function of s alone. Then both r^ k and 
depend only on s, which in turn entails that cM^ fc (s)/<9s, appearing in the 
dual equation of (19), is a function of s alone, as well. However, the phase of 
the exponential in this equation does depend on T. Then it follows from the 
Riemann-Lebesgue lemma [21] that Qnk( s ) § oes to zero in the limit T — > oo. 
This is the argument of Sarandy et al. [15], which allows us to keep only 
the first Pnk( s ) terms under the summation in the dual equation of (17) 
and then uses the adiabatic condition (12) for the a-system, which in turn 
guarantees that the adiabatic approximation is accurate. 
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Next let us focus our attention on the Hamiltonian of the dual b-system 
H b , while we still assume, as above, that H a is T- independent. One can 
see through the relation (14) that rj|^ depends on T on a rapid oscillatory 
manner and through the relation (10) that this applies also to dA b nk (s) / ds. 
Namely, we obtain after some algebra 

<%2 - e "W^ (mjK.) - Kw/&W) j (21) 

where the notation g^ k (s) and tJ^(s) refer to their T- independence. By 
substitution of this equation into definition (19), in the limit T — ► oo we 
obtain 

Q^(*.r)~iT f S rl a k (s')ds' . (22) 

Thus, provided that |r|^| in the above equation is not the constant zero 
function, Q b nk {s,T) may go to infinity in the large T limit. The reason that 
for T — > oo, Q b nk behaves extremely differently than its counterpart Q^ k 
(which rather tends to zero), can be understood by the fact that (contrary 
to dA^h/ds) dA b nk /ds is not a function of s alone, and moreover it depends 
on T in a fast oscillatory manner (see equation (21)). Hence, in the limit 
of T — > oo, dA b nk /ds does not converge for any s £ [0, 1] and in turn the 
Riemann-Lebesgue lemma is not applicable, unlike in the a-system. 

2.5 T-dependence of the b-Hamiltonian 

It is also possible to prove that when H a depends on s only (as we as- 
sumed above), then H b depends also on T, and consequently H b can not be 
written in the only terms of s = t/T. The proof proceeds by reductio ad 
absurdum. Suppose that H b is T- independent. Hence the eigenvalues E\ 
(or equivalently g b nk ) and the eigenstates \E^) of H b are also independent 
on T. According to definition (4) the T-independence also applies to T b k , 

and by virtue of relation (8) it applies to rj]^ as well. However, due to def- 
inition (10) A b nk and then dA b nk /ds should be also T- independent, which is 
the contradiction since dA b nk /ds under (21) depends on T. 

2.6 Extended case 

So far, we have examined the b-system in the large T limit, with the con- 
straint that the Hamiltonian of its dual a-system is independent on T. Now 
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we turn to a more general situation, by imposing the following pair of con- 
ditions on the a-system: 1. it satisfies the adiabatic condition (12), 2. |r|^| 
is strictly greater than zero Vs € [0, 1]. As we can see, now the possible In- 
dependence of the a-system is not ruled out. In the following we prove that 
if these conditions are met, then P^ k is much smaller than Q b nk and then 
P b k can be neglected in the first order approximate solution with respect 
to Qnk- Let us remark that this situation is opposite to the one, discussed 
before in the T-independent a-system, when Q^ k were negligible in relation 
to P^ k . Consequently, let us analyze the fulfillment of 



> 1 . 



(23) 



Subtraction of (17) from (20) and further algebra lead just to the quantity 

iTf*Tl a k (s>,T)ds> 



Q b nk^T) 









1 + 



Pnk(s,T) 



,n / k 



(24) 



The term P b k in (18) can be further written with the aid of (14) and with 
the use of g b nk (s,T) = —g% k (s,T), to obtain 

P b k (s,T) = -(A a nk (0,T)+A a nk (s,T)) . (25) 

Inserting this result back to equation (24) entails the formula 

T 



Qi k {s,T) 




P b nk^T) 





1 



(A« k (0,T)+A« k (s,T)) 



' f S rlt(s',T)ds' 
J 



(26) 



Due to the pair of conditions prescribed for the a-system above, the usual 
adiabatic conditions (12) are satisfied, implying \T/(A® k (0) + A nk (s))\ » 1 

in formula (26) and rj|^(s) is strictly greater then zero Vs G [0,1], hence 
the rhs of (26) is large enough so that inequality (23) holds true and in turn 
P^ k can be neglected in (17). In this case subtracting again (17) from (20) 
yields 

Q b nk (s,T)^ iT f rl a k {s\T)ds' , 



(27) 

which formula can be considered as a T-dependent extension of (22). Thus 
now we have two expressions for Q nk : an exact formula (19) expressed by the 
quantities of the b-system, and an approximate formula (27) defined in terms 
of the non-adiabatic couplings of the a-system, where the approximation 
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is valid provided the a-system satisfies the adiabatic conditions (12) and 
besides t|^(s) is not an extreme small quantity Vs G [0,1]. It is easy 
to construct an a-system which meets the above requirements. However, 
according to the approximate formula (27), and the pair of constraints which 
guarantee the accuracy of (27), Q n k/T is non- vanishing. This implies that in 
the first order solution (17) the amplitudes <j)\ are not constant in time, and 
in turn the evolution of the dual b-system cannot be considered adiabatic. 
Next let us examine, whether the traditional adiabatic conditions, despite of 
the invalidity of the adiabatic approximation, are fulfilled in this b-system. 
According to definition (10), relation (14) and g b nk = —g® k , we obtain 

\A b nk (s,T)\ = \A a nk (s,T)\. (28) 

Thus the b-system satisfies the traditional adiabatic condition (12) only 
when the a-system satisfies it [14] . 

We can summarize that although in the b-system (constructed by rela- 
tion (13) and by imposing the double conditions above) the standard adia- 
batic conditions are satisfied, due to the nonzero values Q b nk /T the state of 
the b-system is not evolving adiabatically. Therefore it represents an exam- 
ple that the standard, traditional adiabatic conditions may not guarantee 
adiabaticity (Tong et al. [14]). On the other hand, we argued that the 
Hamiltonian of this b-system can not be written in terms of t/T quanti- 
ties, unlike the case when the standard adiabatic conditions are sufficient in 
guaranteeing adiabaticity [15, 19]. 



3 Illustration of the Theory 

In the following the general analysis of the previous section is supported by 
a simple example to demonstrate that the non- vanishing Q b nk /T terms has 
a key role, in that the standard adiabatic conditions are not sufficient to 
guarantee the validity of the adiabatic theorem. To this end imagine that 
the a-system is a spin-half particle at rest at the origin in the presence of a 
constant magnitude magnetic field rotating in a plane at constant angular 
velocity to. This type of system was also considered by refs. [22, 23]. The 
corresponding effective Hamiltonian in the rest frame is 
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where ujq is defined by the magnetic moment of the spin and the strength of 
the magnetic field. The corresponding eigenvalues are Ef(s, T) = —E!j;(s, T) = 
ujq/2. Identifying T = 2-k/uj and switching to normalized time s by the 
transformation t = sT, we obtain for the normalized eigenspinors 

\Ei(s,T) = ^ ( e ~™ s ) , \E%(s,T) = ^ ( ) , (30) 

so that neither the Hamiltonian H a nor the corresponding eigenvectors are 
T-dependent. The local phases of the eigenstates has been chosen so that 
they obey the parallel transport law (r^„(s,T) = for n = 1,2). In this 
case by direct substitution of (30) into definition (4) tJ" = rfi = —in. 
Calculating the Hamiltonian H h of the dual system from (13) to the first 
order in the small quantity ^ one finds 

uj I cos looTs i sin ujqTs \ . . 
2 I — i sin ujqTs — cos ujoTs J 

Evidently H b has an irremovable T-dependence. Let us now treat the b- 
system perturbatively, by the general method given above. We consider 
that the dual b-system is initially in its instantaneous \El(0,T)) eigenstate, 
and calculate the relevant Q2i(s,T) term both exactly and approximately 
by the mean of equations (19) and (27), respectively. The approximate solu- 
tion according to the previous section is good if the adiabatic condition (12) 
holds for the a-system and that |r|"| is nonzero. In our two- level model the 

absolute value of tJ" is tt, whereas condition (12) implies uj <C ujq. If T is 
big enough so that uj <C ujq fulfills, then the exact solution is approximated 
sufficiently well. However, in the concrete calculations both the exact and 
approximate solutions yield Q2i{s,T)/T = —its, which surprising result is 
due to the fact that the above pair of conditions are in fact not necessary for 
the approximation would be faithful (i.e. they are merely sufficient condi- 
tions). Thus the formula Q\\jT = —its tells us that regardless of the value 
of T, Q21/T may become relevant provided s is greater then zero. 

Let us now calculate how faithful the adiabatic approximation is in this 
case. The adiabatic approximation is acceptable in the interval t £ [0, T] 
if the overlap \^(t)\Ei(t))\ 2 = |(V>(s,T)|£i(s,T))| 2 ~ 1 for s £ [0,1] pro- 
vided the state of the system is prepared initially in the eigenstate \Ei(0)). 
However, 

IWs,T)|£i(s,T))| 2 = 1 - ms,T)\E 2 (s,T))\ 2 = 1 - |0 2 ( S ,T)| 2 , (32) 
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where the first equality is valid for two-level systems, and the expansion for- 
mula (5) was used to obtain the second equality. For our concrete example, 
according to (20) in first order approximation the amplitude 4> h 2 reads 

4(s,T) = - f S r^(s',T)ds' = ins , (33) 
Jo 

and by returning to the original time variable t = ^ we can write this 
further as 

|(^(t)|^(t))| 2 = l-(y) 2 • (34) 

In our simple example rj" does not depend on s and equation (16) can be 
easily solved (without applying perturbation theory) to obtain the exact 
solution 02(s,T) = sin(7rs) = sin(^), which yields the exact fidelity 

|^)|^(t)| 2 = l-sin 2 ^ (35) 

for our particular example. If we compare (35) with (34), one can see that 
the approximate formula (34) reproduces the first term of the Taylor series 
expansion of the sine function in the exact formula (35). Performing higher 
order perturbation by means of an iterative solution of equation (15), the 
resulting series will approach more and more accurately to the exact expres- 
sion for the overlap. This implies that in order to solve (15) by perturbation 
method one has to take into account at each iterative step non-negligible 
off-diagonal terms (where in the first order the corresponding term is Q h nk 
for n / k). These terms when we go beyond the first order approximation 
accumulate and according to the exact formula (35) for the overlap (alterna- 
tively see Tong et al. [14]) for sufficiently large t results in a great deviation 
from the ideal overlap 1, pertaining to the pure adiabatic evolution. 

In summary, this simple example shows that correct to first order it 
is indeed Q n k , (n = 1, k = 2) term in (17) that is responsible for the 
breakdown of the adiabatic approximation, despite the fulfillment of the 
traditional adiabatic conditions. 



11 



4 Test of the novel adiabatic condition 



Finally, we test a new quantitative adiabatic condition, proposed by Ye et 
al. [16]. The new quantity, introduced in ref. [16], looks as 



A nk (s,T) 



g n k(s,T) 



1 d,argTH k (s,T) 
T ds 



(36) 



and the new quantitative adiabatic condition is encoded in the requirement 
[16] 

\A nk {s,T)\ <T, k + n, s€[0,l], (37) 
in order the adiabatic approximation be faithful in the time range t 6 [0, T\. 



By substitution of the relations (14) and g° nk = —g^ k 
systems, the following simpler condition is given for the b-system, 



between the dual 



A b nk (s,T) 




ll*>/ 
T nk( S 


T) 




r JIfc( S > T ) 


T 




Tg b nk (s,T)- 


dargr^T 
ds 




9ax gT >( S ,T) 
ds 



< 1 . (38) 



Next the above condition is applied for our special two-level b-system, dis- 
cussed in the previous section. Since rj" (s,T) = — iir, its argument remains 
constant in time, and therefore quantity \A b nk (s, T)/T\ is infinite for any 
s G [0, 1] and in turn (38) is not fulfilled. This implies that the novel 
adiabatic condition of Ye et al. [16] indicates correctly that the adiabatic 
approximation is violated in the b-system, whereas the traditional adiabatic 
condition fails to detect it. 



5 Conclusion 

In conclusion we investigated through a first order perturbative solution the 
state evolution of a dual pair of systems, which construction was considered 
first by Marzlin and Sanders [13]. The a-system is assumed to satisfy the 
usual adiabatic criteria, and the b-system performs a reversed time evolu- 
tion with respect to its dual pair. In this case we argue that in the b-system 
due to the fast oscillations of the off-diagonal non-adiabatic coupling terms 
(t^), the terms Q b nk defined by (19) may become large and dominate in the 
first order solution of (15). This results in the breakdown of the adiabatic 
approximation in the b-system, although the standard adiabatic conditions 
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are satisfied in it. Thus, we could find via perturbative treatment the un- 
derlying reason why the traditional adiabatic conditions are insufficient in 
the b-system. It is further argued that in this case the Hamiltonian of the 
b-system cannot be written in the form of H(t/T). The above findings are 
illustrated on a simple two-level model as well, whereon a new adiabatic 
condition was also tested and showed that in this particular case it performs 
well. 
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